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Abstract 

We describe the C program mori.x. It is part of PALP, a package for analyzing lattice 
polytopes. Its main purpose is the construction and analysis of three-dimensional smooth 

■ Calabi-Yau hypersurfaces in toric varieties. The ambient toric varieties are given in terms 
£S| of fans over the facets of reflexive lattice polytopes. The program performs crepant star 

■ triangulations of reflexive polytopes and determines the Mori cones of the resulting toric 
varieties. Furthermore, it computes the intersection rings and characteristic classes of 

\Q ■ hypersurfaces. 
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1 Introduction 

PALP PQ, a package for analyzing lattice polytopes, has the ability to construct, manipulate 
and analyze lattice polytopes. The present program mori.x, which is also included in the new 
version of PALP (starting from release 2.0) available at [2J, adds further functionalities con- 
cerning Calabi-Yau (CY) three-fold hyper surf aces. The smoothness of the CY hypersurface 
is achieved by triangulating the reflexive polytope, i.e. appropriately resolving the ambient 
space. This functionality is optimized for the CY hypersurface case. Generic CY hypersur- 
faces avoid point -like singularities of the ambient space as well as divisors that correspond 
to interior points of facets. Consequently, our algorithm performs crepant star triangulations 
only up to such interior points Using these triangulations, the SR-ideals and the inter- 
section rings are computed. The latter are determined with the help of SINGULAR [5], a 
computer algebra system for polynomial computations. Furthermore, we have implemented 
the Oda-Park algorithm [51[7j to find the Mori cone of the ambient space. 

The program can also analyze arbitrary three-dimensional hypersurfaces embedded in 
the ambient toric varieties. It is capable of computing the intersection ring, characteristic 
classes and Hodge numbers. As we do not triangulate the polyhedron completely, non-CY 
hypersurfaces can also contain singularities. Hence, the intersection numbers of non-CY 
hypersurfaces determined by mori.x are not integers in general. 

The triangulations, which are at the heart of this program, are subject to two types of lim- 
itations. First, mori.x performs crepant star triangulation of the polytope by subdividing its 
facets into simplices. This procedure is, however, only implemented for those non-simplicial 
cones which give rise to at most three-dimensional secondary fans. Second, the limitations of 
PALP apply, as mori.x is a part of this package. PALP works with fixed precision; maximal 

1 Complete triangulations of arbitrary polytopes can be performed with the program TOPCOM [3], which 
is also included in the open-source mathematics software system SAGE [4]. SAGE also contains various tools 
for handling toric varieties. 
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dimensions, maximal number of polytope points, vertices and facets need to be fixed before 
compiling. By editing the file Global. h, the user can change the default values if needed. See 
the PALP documentation pQ for further details. 

This note is meant as a manual for mori.x, an introduction to the mathematical back- 
ground is beyond its scope. We refer the reader to [8,9j for introductions to toric geometry 
and the construction of CY hypersurfaces using polytopes. The structure of the paper is as 
follows. In section 2 we explain the I/O formats, present the options of mori.x and describe 
the basic functions of the main program. This section should be used as a quick reference 
for the various functionalities. Aspects which deserve further explanations are presented in 
section 3. There, we also discuss the resolution of singular compact toric spaces via (crepant 
star) triangulations. Section 4 describes the structure of the program. The principal routines 
are listed and commented on in the header file Mori.h. 

Acknowledgments: We dedicate this work to Maximilian Kreuzer who passed away be- 
fore it was completed. He contributed the core routines to the source code of mori.x. The 
present publication is an occasion to savor his work. We hope that we have been able to 
realize the original spirit in which he started this project. 

We would like to thank Harald Skarke for valuable advice in the completion phase of 
the program, and Christoph Mayrhofer for his comments and contributions to the source 
code at an early stage. Furthermore, we would like to thank Johanna Knapp, Volker Braun, 
Maximilian Attems, Ching-Ming Chen and Michal Michalcik for discussion and comments. 

The work of N.-O. Walliser was supported by the FWF under grant P21239 and 1192. 
The work of A. P. Braun was supported by the FWF under grant P21239 and P22000. 

2 I/O formats and options of mori.x 

In this section we explain the I/O formats of mori.x and give an overview of its options. This 
program is part of the latest version of the package PALP that is available at the web page [2 J 
under the GNU license terms. There, the reader can find a detailed installation description. 
In short, type "make all" in the command line to compile the entire suite of programs. 
Otherwise, type "make mori" to compile mori.x only. Consult pQ for the documentation of 
the other functions of the package: poly.x, cws.x, class. x and nef.x. 

Let us start with the help screen. It provides essential information about all the function- 
alities of the program (letters in boldface denote the command line input): 

$ mori.x -h 

This is ''mori.x'': performing crepant star triangulations of a polytope P* in the N-lattice 
computing the Mori cone of the corresponding toric ambient spaces 
determining the intersection rings of embedded (CY) hypersurfaces 
Usage: mori.x [-<Option-string>] [in-file [out-file]] 
Options (concatenate any number of them into <Option-string>) : 
-h print this information 

-f use as filter 
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-g general output: triangulation and Stanley-Reisner ideal 

-I incidence information of the facets (ignoring interior points of facets) 

-m Mori generators of the ambient space 

-P IP-simplices among points of P* (ignoring interior points of facets) 

-K points of P* in Kreuzer polynomial form 

-b Hodge numbers and Euler number 

or arithmetic genera and Euler number if combined with -H 
-i intersection ring 

-c Chern classes of the (CY) hypersurface 

-t triple intersection numbers 

-d topological information on toric divisors and del Pezzo conditions 

-a all of the above except h, f, I and K 

-D lattice polytope points of P* as input (default CWS) 

-H arbitrary hypersurface class 'H = cl*Dl + c2*D2 + . . . ' as input (default CY) 

Input: coefficients 'cl c2 
Input: 1) standard input format: degrees and weights 'dl wll wl2 . . . d2 w21 w22 
2) alternative (use -D) : 'd np' or 'np d' (d=Dimension, np=#[points]) 

and (after newline) np*d coordinates 
Output : as specified by options 



As the program works with reflexive polytopes only, the input of the program must specify 
such a polytope. There are two input formats available: one either gives a combined weight 
system (CWS) or directly provides a matrix of lattice points whose convex hull is the polytope. 
The input is checked for reflexivity, see [HQS] f° r details on how a CWS is converted into a 
polytope. 

The construction of a CY hypersurface uses a dual pair of reflexive lattice polytopes. One 
polytope determines the CY hypersurface, whereas its dual specifies the ambient toric variety. 
In this work, we are mainly interested in the iV— lattice polytope which gives rise to the fan of 
the ambient toric variety. We follow the conventions of the literature and PALP in referring 
to this polytope as the dual polytope P* . 

Let us consider a P 1 fibered over a P 3 as an example: 

$ mori.x -P 

Degrees and weights 'dl wll wl2 . . . d2 w21 w22 ...': 

5111110 2 8 8 1 1 

4 7 points of P* and IP-simplices 

1 8 8 -1 

8 18 8-100 

8 8 18-100 

8 8 8 1 1 -1 
#IP-simp=2 

1 1 1 8 1 1 5=d codim=8 

8 8 8 18 1 2=d codim=3 

Note that the ordering of the CWS input is not obeyed by the output of lattice points. Once 
the order is displayed, however, it is fixed and determines the labeling of toric divisors in any 
further output. 

An alternative way to provide the input is to type lattice polytope points directly. In this 
case, one has to use the parameter "-D" . Let us reconsider the example above: 
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$ mori.x -DP 

'#lines #colums' (= 'PolyDim #Points' or '#Points PolyDim'): 
4 7 

Type the 28 coordinates as dim=4 lines with #pts=7 columns: 
19 0-1 000 
10-1 000 
1-1 
1-110 

4 7 points of P* and IP-simplices 



1 








-1 














1 





-1 

















1 


-1 




















1 


-1 


1 






#IP-simp=2 

1 1 1 1 1 5=d codim=0 
11 2=d codim=3 



The order of the lattice points displayed is the same as the order set in the input. This type 
of input should be preferred if one wants to control the order of the polytope points. Note 
that the IP-simplices among points of an M-lattice polytope P can be obtained from the 
lattice points of P* with the option "-P" of poly.x. Since we focus on the iV-lattice polytope 
P*, this option is suppressed in mori.x for the sake of simplicity. 

The rest of this section contains a detailed description of the options listed in the help 
screen. If no flag is specified, the program starts with the parameter "-g". By default, the 
program considers a CY hypersurface embedded in the ambient toric variety. The option 
"-H" has to be used in order to consider non-CY hyper surf aces. 

-h This option prints the help screen. 

-f This parameter suppresses the prompt of the command line. This is useful if one wants 
to build pipelines. 

-g First, the triangulation data of the facets is displayed. The number of triangulated 
simplices is followed by the incidence structure of the simplices. The incidence infor- 
mation for each simplex is encoded in terms of a binary number: there is a digit for 
each polytope point; a 1 denotes that the point belongs to the simplex. Second, the 
Stanley-Reisner ideal is displayed: the number of elements of the ideal is followed by 
its elements. Each element is denoted by a binary number as above. 

-I The incidence structure of the facets of the polytope P* is displayed. Interior points of 
the facets are neglected. 

-m The Mori cone generators of the ambient space are displayed in the form of a matrix. 
Each row corresponds to a generator. The entries of each row are the intersections of the 
generator with the toric divisor classes. The Oda-Park algorithm is used to compute 
the generators. 
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-P Those IP-simplices whose vertices are lattice points of the polytope, but not inner points 
of the facets, are displayed. 

-K The Kreuzer polynomial is displayed, see section [3721 for details. The number of points 
in the interior of facets is shown as intpts. The multiplicities of the toric divisors are 
displayed as multd if they are greater than one. Furthermore, the Picard number of the 
CY hypersurface is computed and printed as Pic. 

-b We distinguish two cases. For the Calabi-Yau hypersurface (default) case, the Hodge 
numbers h ' , h 2,1 and the Euler characteristic are displayed. If the input is an arbitrary 
hypersurface (see option "-H"), the zeroth and first arithmetic genera are displayed 
instead of the Hodge numbers. As a non-CY hypersurface can be singular, care is 
needed for an interpretation of the results; see section 13.11 for more details. 

-i This option displays the intersection polynomial restricted to a CY hypersurface. The 
polynomial is displayed in terms of an integral basis of the toric divisors. The coefficients 
of the monomials are the triple intersection numbers in this basis. This option can also 
be used together with "-H" to perform this task for non-CY hypersurfaces. 

-c The Chern classes of the hypersurface (CY or non-CY) are displayed in terms of an 
integral basis of the toric divisors. 

-t The triple intersection numbers of the toric divisors are displayed. When the intersection 
numbers are equal to zero, they are not shown. 

-d This option displays topological data of the toric divisors restricted to the (CY or 
non-CY) hypersurface. The Euler characteristics of the toric divisor classes and their 
arithmetic genera are shown. The toric divisor classes are tested against necessary 
conditions for del Pezzo surfaces. The following data is listed: the del Pezzo candidates 
preceded by their number (their type is given in brackets) and those del Pezzo candidates 
that do not intersect other del Pezzos. 

-a This is a shortcut for "-gmPbictd". 

-D This tells the program to expect a matrix of lattice points of the polytope P* as the 
input. This is useful if one wants to set and control the order of the toric divisors. 

-H Using this option, one can specify a (non-CY) hypersurface. The user determines the 
hypersurface divisor class H = CjZ?j in terms of the toric divisor classes Di by typing 
its coefficients q . The hypersurface can then be analyzed by combining "-H" with other 
options, as described above. Just using "-H", the program runs "-Hb". See section [3711 
for an example. 

Note that the options "-b, -i, -c, -t, -d, -a, -H" need SINGULAR to be installed. 
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3 Supplementary details 



We designed this program with an eye for applications to string theory model building, in 
particular the analysis of CY three-folds. In this context, our motivation was to implement the 
construction of smooth CY hypersurfaces embedded in toric varieties, using the construction 
of Batyrev [11] . Here, the starting point is a dual pair of reflexive polytopes which determines 
both the ambient toric variety and a CY hypersurface. 

3.1 Triangulations and point— like singularities 

Batyrev [11] has shown that any four-dimensional reflexive polytope gives rise to a smooth 
CY hypersurface after triangulation. The ambient toric space is not necessarily smooth, 
as some of the cones might have a volume greater than one even after triangulation. This 
leads to point-like singularities, which, however, do not meet a generic CY hypersurface. 
The reflexivity of the polytopes sensibly simplifies the triangulation procedure. The induced 
simplices of lower dimension do not contribute any further singularities. 

mori.x performs crepant star triangulations of a polyhedron ignoring points in the interior 
of facets. Even though this introduces further point-like singularities into the ambient toric 
variety, these are also avoided by generic CY hypersurfaces. 

Polytopes can be triangulated by subdividing the secondary fans of its non-simplicial 
facets [12] I13j . This triangulation algorithm is implemented in mori.x for up to three- 
dimensional secondary fans. The program exits with a warning message if the subdivision is 
not properly completed. 

Consider the following CWS: 

$ echo -e '8 411118 6 3 18 18 1' | mori.x -fPI 

4 8 points of P* and IP-simplices 

3 1 8 8 8 -1 1 
-1 1 
3 1 -1 1 
-4 1 1 -1 

#IP-simp=2 

1110 14 8=d codim=0 
10 110 3 6=d codim=l 
Incidence: 101011 001111 111110 110101 011101 111001 100111 

This system describes a four-dimensional lattice polytope with eight points, six of which are 
vertices. We label the column vectors with v±, . . . ,v$. Point vj lies in the interior of the third 
facet. The incidence data show the intersections of the six polytope points with the seven 
facets. The third facet contains the five points v\, . . . ,v$, hence it is not simplicial and we 
have to triangulate it. 

$ echo -e '8411118 6318181' | mori.x -fg 

9 Triangulation 

101011 001111 110101 011101 111001 100111 011110 110110 111010 
2 SR-ideal 
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101100 010011 
8 Triangulation 

101011 001111 110101 011101 111001 100111 101110 111100 
2 SR-ideal 
010010 101101 

The program performs the two possible crepant triangulations of the facet (12345). The first 
result yields the three simplices (13425), whereas the second gives the two simplices (12534) 
(in this notation the hat indicates that one of the points is dropped). Nevertheless, the two 
resolutions give the same CY intersection polynomial: 

$ echo -e '8411118 6318181' | mori.x -fi 

SINGULAR -> divisor classes (integral basis Jl ... J2): 
dl=Jl, d2=-3*Jl+J2, d3=Jl, d4=4*Jl-J2, d5=-3*Jl+J2, d6=J2 
SINGULAR -> intersection polynomial : 
2*J1"3+188*J2"3+8*J1"2*J2+38*J2"2*J1 

SINGULAR -> divisor classes (integral basis Jl ... J2): 
dl=Jl, d2=-3*Jl+J2, d3=Jl, d4=4*Jl-J2, d5=-3*Jl+J2, d6=J2 
SINGULAR -> intersection polynomial: 
2 * J 1" 3 + 1 8 * J 2 - 3 + 8 * J 1" 2 * J 2 + 3 * J 2 " 2 * J 1 

dl,..., d6 denote the toric divisors corresponding to the lattice points vi,...,v$. There 
are two independent divisor classes. Indeed, mori.x expresses the intersection polynomial in 
terms of the integral basis J\ = D\ = D% and J2 = Dq. 

Let us take a closer look at non-CY hypersurfaces. The reader is warned: for these 
cases smoothness is not guaranteed anymore, so that the intersection numbers can become 
fractional. Some choices of the hypersurface equation may intersect point-like singularities 
not resolved by the triangulation. Consider e.g. the hypersurface divisor class H = D\ + Dq. 
Remember that the order in which mori.x expects the coefficients of the hypersurface divisor 
class is fixed by the polytope matrix and not by the CWS input. Hence, the correct input for 
H is the string 1 © © 1. 

$ mori.x -H 

Degrees and weights 'dl wll wl2 . . . d2 w21 w22 
8411110 6310101 

Type the 6 (integer) entries for the hypersurface class: 
1 1 

Hypersurface degrees: (5 4 ) 
Hypersurface class: l*dl l*d6 

SINGULAR -> Arithmetic genera and Euler number of H: 
chi_0: 35/32 , chi_l: 143/32 [ -27/4 ] 
SINGULAR -> Arithmetic genera and Euler number of H: 
chi_0: 29/27 , chi_l: 128/27 [ -22/3 ] 

To calculate these quantities, the program determines the characteristic classes of the divisors 
using adjunction. It then performs the appropriate integration with the help of the triple 
intersection numbers. The fractional results of the arithmetic genera and the Euler number 
in our example indicate that the intersection polynomial has fractional entries. This happens 
because the singularity of the ambient toric variety descends to the hypersurface H. In the 
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first triangulation, the simplex (1235) still has volume four. In the second triangulation, 
there are two simplices with volume three: (1234) and (1345). Indeed, any hypersurface in 
the divisor class H = D% + Dq is forced to pass through the corresponding singularities. 

3.2 The Kreuzer polynomial 

The Kreuzer polynomial encodes lattice polytope points in a compact form. The number of 
variables equals the dimension of the polytope. Each lattice point gives rise to a Laurent 
monomial in which the exponents of the variables are the coordinates. Vertices and non- 
vertices are distinguished by coefficients "+" and "-" respectively. Points in the interior of 
facets are ignored. Consider the example presented above. 

$ echo -e '8411118 6 3 16 18 1' | ./mori.x -fK 

KreuzerPoly=t_l"3t_3"3/(t_2t_4"4)+t_l+t_2+t_4+t_3+t_4/(t_lt_3) ; intpts=l; Pic=2 

Note that negative coordinates are always displayed by putting the variables in the denomi- 
nator. 

3.3 The Mori cone of the ambient space 

For toric varieties, mori.x uses the algorithm of Oda and Park [6K7] to compute the Mori cone. 
The generators of the Mori cone are given in terms of their intersections with the toric divisors. 
For singular toric varieties, the Picard group of Cartier divisors is a non-trivial subgroup of 
the Chow group, which contains the Weil divisors. Hence one can consider the Kahler cone, 
which is dual to the Mori cone, as a cone in the vector space spanned by the elements of 
either the Picard or the Chow group. The program mori.x only deals with simplicial toric 
varieties, for which the Picard group is always a finite index subgroup of the Chow group [8]. 
Hence the Cartier divisors are integer multiples of the Weil divisors and this ambiguity does 
not arise. 

3.4 The topological data of toric divisors 

Using option "-d", mori.x computes the arithmetic genera of the toric divisors restricted 
to the embedded hypersurface and determines the del Pezzo candidates among them. The 
program checks the del Pezzo property against two necessary conditions: First, for a del Pezzo 
divisor S of type n, the following equations should hold: 

[ Cl (S) 2 = 9-n, [c 2 (S)=n + 3 => Xo (S) = [ Td(S) = 1 . (1) 
Js Js Js 

Here, Td(S) denotes the Todd class of S, which gives the zeroth arithmetic genus of S upon 
integration. This test also allows to determine the type of the del Pezzo surface in question. 

The second necessary condition comes from the fact that a del Pezzo surface is a two- 
dimensional Fano manifold. Hence, the first Chern class of S integrated over all curves on S 
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has to be positive: 

D i nS , nci(5) >0 VDi-.Dij^S, (2) 

This condition would be sufficient if we were able to access all curves of the hypersurface. 
In our construction, however, we can only check for curves induced by toric divisors. This 
functionality was added to carry out the analysis of base manifolds for elliptic fibrations 
in HI]. 

4 Structure of the program 

mori.x is part of the new releases of PALP (starting from version 2.0). The general structure 
of the package has not been changed, only some new files have been added. Hence, all general 
annotations to the package in pQ remain valid. In this section, we provide an overview of the 
composition of mori.x and discuss its dependencies on pre-existing files. 

The source code of mori.x is contained in the program files mori.c, MoriCone.c, Singu- 
larlnput.c, and the header file Mori.h. Makefile reflects the dependencies of the program. 
The compilation supports 32 as well as 64 bit architectures. One can adjust the compilation 
parameters according to ones needs. The optimization level is set at -03 by default. 

mori.c contains the main and the help information routines. Further, basic manipulations 
(completion, calculation of facet equations,...) of the polyhedron are performed with the help 
of core routines from Vertex. c. In particular, reflexivity of the polytope is checked. 

MoriCone.c is at the heart of mori.x. After determining the non-simplicial facets, their 
triangulation is performed by the routine GKZSubdivide. This function identifies the maximal 
dimensional secondary fans of the facets and makes a case-by-case triangulation depending 
on their dimensions. This function is only implemented for secondary fans up to dimension 
threej§ Once the subdivision is accomplished, the program determines the Stanley-Reisner 
ideal (StanleyReisner) and computes the Mori cone (Print_Mori). Furthermore, it finds a 
basis of the toric divisor classes. 

Singular Input, c is the interface to SINGULAR [5] . The latter is a very efficient computer 
algebra system for computations with polynomial rings. In Singular -Input, c, the Chow ring 
is determined from the Stanley-Reisner ideal, the linear relations among the toric divisors, 
and a basis of the toric divisors of the triangulated polytope. This data is put together by 
HyperSurfSingular and redirected to SINGULAR, which then determines the intersection 
ring restricted to the hypersurface and computes its characteristic classes. |§ 

The most important routines of mori.x are documented in the file Mori.h. This header 
file provides a more detailed description of the structure of the program. 

2 The pre-compiler command TRACE_TRI ANGULATION in MoriCone.c enables diagnostic information 
about the triangulation. This data might be of use for the motivated programmer who wants to extend the 
subdivision algorithm. 

3 The input for SINGULAR can be displayed in the standard output of mori.x by turning on the pre- 
compiler definition TEST_PRINT_SINGULARJO in Singularlnput.c. 
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